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MULTIPLE MIDCOURSE MANEUVERS IN 
INTERPLANETARY GUIDANCE 

by 

Abolghassem Ghaffari 


SUMMARY 

An attempt has been made to investigate the number and location of mid- 
course maneuvers in an interplanetary mission. 

The case of a fly-by Jovian mission with a spin stabilized spacecraft is con- 
sidered. The mission objectives are a fly-by within 10 6 km (6 diameters) of the 
planet Jupiter and a flight time from earth to Jupiter of 500 days. 

The spacing-ratio method has been applied and it is shown that a constant 
spacing-ratio of 2.5 for the velocity-sensitivity leads to a near optimum spac- 
ing of correction times. 
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MULTIPLE MIDCOURSE MANEUVERS IN 
INTERPLANETARY GUIDANCE 


INTRODUCTION 

A spacecraft moving in free flight toward some target region in space is 
guided to its final destination by applying one or more small velocity impulse 
corrections at certain times along the trajectory to null the predicted target 
error. The estimate of the target error is achieved by an orbit determination 
process; the required corrections are computed using linear perturbation theory, 
and the impulse is delivered by a rocket motor, which applies an acceleration to 
the spacecraft for a relatively short period of time. The selection of times for 
performing the velocity corrections to the actual trajectory, and the determina- 
tion of what fraction of the predicted target errors is to be nulled by each ma- 
neuver is termed by guidance law or guidance policy. 

The technique for applying a single-impulse velocity correction to the tra- 
jectory of a spacecraft has been developed at Jet Propulsion Laboratory and is 
now almost classic. This technique has led to the Mariner II fly-by mission to 
the planet Venus and also to the Mariner IV fly-by mission to the planet Mars. 

In this case, a suitable single-maneuver time is chosen from preflight studies 
of orbit determination and execution error statistics, and the correction capa- 
bility to be carried aboard the spacecraft is determinated by mapping the covar- 
iance matrix of injection guidance errors to the selected maneuver point [1] to 
obtain the covariance matrix of ’ 've locity-to -be -gained ' ' components. 

The consideration and realization of certain flights with more rigid mission 
objectives and higher accuracies such as Apollo mission, satellite injection 
about planets and soft landing, led us to believe that a single maneuver is inade- 
quate to accomplish the mission objectives with the required degree of accuracy. 

The accuracy that can be obtained from multiple midcourse corrections is 
a considerable improvement over the single maneuver. On the other hand, the 
situation becomes much more complex when more than one maneuver is con- 
sidered, because the future guidance and tracking policy must be considered for 
the application of a maneuver at any given time. 

The target error criterion and the observed accuracy must be defined, as 
well as the bound on the total velocity correction that can be applied. This 
question has been investigated, as an optimization problem by Battin [2], Break- 
well [3], Striebel and Breakwell [4] , Lawden [5], Pfeiffer [1], and others. 
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In general, there would be a sequence of times at which impulsive correc- 
tions might be made. Assuming that the task of the guidance system is to null 
the predicted target error, and that orbit determination data to be taken are 
specified, the sequence of times for performing the corrections is associated 
with the guidance law. The selection of these times is an example of an opti- 
mal stochastic control problem. 

To arrive at criteria for the selection of these times, several analytical 
methods have been proposed. The main techniques developed are the three 
following schemes for the determination of the guidance law. 

1. The variance ratio method, suggested by Battin [2] , in which a criterion 
for the timing of corrective impulses has been proposed. This criterion 
is based on the ratio of the required velocity correction to the uncer- 
tainty in estimated miss distance, assuming that each maneuver nulls 
out the estimated miss distance error. This criterion is not directly 
related to a concept of minimum fuel expenditure. 

2. The spacing ratio method, developed by Breakwell [3], in which one 
seeks the timing of corrective impulses to minimize the average total 
maneuvers under the assumption that: 

a. Each maneuver nulls out the estimated miss distance error. 

b. The error in estimating miss distance is due entirely to an error 
in estimating the instantaneous velocity vector. The optimum 
choice for the spacing of corrections is one which minimizes 

the total velocity correction or, equivalently, minimizes the total 
required fuel. 

3. The minimum error method, proposed by Pfeiffer [1] , in which the 
problem is approached from a dynamic programming point of view in 
order to formulate an adaptive policy that seeks to minimize the mean- 
squared target error, subject to a propellant constraint. In this 
method the maneuver times t . are designated, and the possibility 

of performing a maneuver at each of these times is allowed. 

The purpose of this paper is to investigate the number and location of the 
correction points on a fly-by Jovian mission. The choice of the spacing- ratio 
method is mainly due to its characteristic features. 

The symbols and nomenclature used in this paper are taken from reference 
[6] which are commonly used by the Mission Analysis Office. The nominal tra- 
jectory used for this Galactic Probe is adopted from [7] with the same elements 
and characteristics. 
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STATEMENT AND ANALYSIS OF THE PROBLEM 

Suppose a spacecraft moves in free flight on its way from earth to the planet 
Jupiter, so that its trajectory would pass within 10 6 km of the planet Jupiter. The 
spacecraft attitude will be restricted during its motion, and more specifically 
its attitude control law will be assumed that its spin axis A° (unit vector) iner- 
tially fixed as parallel to the injection velocity AV, Figure 1, (reference 8), that 
is to say: 



( 1 ) 


where 


AV° 


AV 

I AV | 


In fact, if proper precautions are taken to prevent decay of spin rate, spin sta- 
bilization has the most general application of all passive systems. 



0 = COS - 1 (A° • r° ) SPIN AXIS - EARTH ANGLE 
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Figure 1 -Jupiter Transfer Geometry 
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Midcourse maneuvers are also restricted to be applied along the spin axis 
A°. 


Except in the immediate vicinity of the spheres of influence of the earth and 
the planet Jupiter, the nominal trajectory of the spacecraft is essentially a heli- 
ocentric ellipse. This problem is three-dimensional, but we assume the orbits 
of earth and Jupiter and transfer ellipse are coplanar, and the treatment is two- 
dimensional. 

The actual transfer trajectory is somewhat different from the nominal tra- 
jectory due to injection errors and therefore the spacecraft would miss the de- 
sired terminal target by a miss distance B(t n j) if no corrective action were 
taken. 

It is desired to apply at time t n a corrective velocity impulse AV n at junc- 
tion point C n , Figure 2, to reduce the error at the terminal point. 



Figure 2-Interplanetary Trajectory 
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By applying a linearized error theory one finds [3] 1 the expression of the 
computed corrective velocity impulse AV n as follows: 

(^r) • (“-) +i ^> = °- < ! 


where B(t n _j) denotes an estimate of the miss distance B(t n _ 1 ) prior to the 
nth impulse and subsequent to the (n - l) th impulse. The partial derivative is to 
be evaluated along the nominal trajectory on which there are no impulses. 

It is convenient to use the subscript n and set, in the sequel, 


B(t ) = B . 

v n 7 n 

The miss distance B n _ 1 is not usually correctly estimated, that is to say, 
a correction at the junction point C n still leaves us with a miss distance B n 
which may be reduced by further corrections. If the estimated B n-1 differs 
from the actual B n _j by an amount of e n-1 , then 


B . = B , - e , 

n - 1 n-1 n - 1 


( 3 ) 


If y° and A V° denote the unit vectors along the vectors Bb/BV and AV n respec- 
tively, then the equation (2) can be written: 



^ 

BB 

AV n 

BV 



cos (7 0 


AV n °) + 


B n - X =0. 


( 4 ) 


The minimum of AV n requires the condition 


cos (7° , Av n °) = 1 


( 5 ) 


i.e., the velocity-sensitivity vector BB/BV and the veloc ity correction vector 
Av" are parallel. They have the same direction if B n _ x < 0 and opposite di- 
rection if B n _ i > 0 according as the passage of the spacecraft is at right or 
left of Jupiter. The computed velocity correction is therefore 


iFor a detailed and complete analysis of the spacing-ratio method, see ref. 


[3]. 
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AV 


Bl 

BV 


( 6 ) 


Conditions (1) and (5) leads us to conclude the following optimal constant- 
attitude criterion: 


The constant -attitude thrust optimization leads to constant-attitude velocity- 
sensitivity . 

If the magnitude of the actual mechanized velocity correction is 


Av » = Av „ + v; , (7) 

where V' is the velocity mechanized error, then the expression of the miss 
distance after nth correction becomes 


BB_ 


B = B , + (AV + V') — IL . 

n n-1 v n n 7 ~\y 


( 8 ) 


Taking into account of the equations (2) and (3), the condition (8) becomes 


B = e„ 

n B, 


n - 1 


bb 

+ V' — 
n BV 


(9) 


From equations (6), (7), and (9) we obtain 


B 


AV = -- n -~- 1 . - + V' . 
n BB n 

n 

"bv 


( 10 ) 


Replacing n by n - 1 in equation (9) and using equations (3) and (10) we now 
obtain 


AV = V* + 

n n 


BB 


BV 


-l 


v,*v. 


3B n 1 

_ v' n l 

n-1 3 V 


(ID 


Equation (11) gives us the expression for the nth velocity correction and 
shows its connection to errors at and prior to the nth junction point. 
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Now the main concern is, given t N and an RMS launching error cr Bo , choose 
the correction or junction points ,C 2 ,. . ..C^, Figure 2, as to achieve the re- 
quired mission objectives with a minimum total velocity correction and hence a 
minimum fuel expenditure for corrective thrusts. In other words, given t N and 
the standard deviation a B , choose a sequence of times {t ; n = 1, 2, . . . ,n} , 
with the integer N not specified such that the sum of the magnitude of the veloc- 
ity corrections 


N 

F I AV„| (12) 

n=l 

or more workable and closely related expression 

N 

S*=2]E[|AV n |] (13) 

ns 1 

is as small as possible. N denotes the total number of corrective thrusts, and 
| A V n | is the magnitude of the impulse applied at time t n . The symbol E [ | A V | ] 
indicates the mean value or expected value of the random variable | AV n |. 

If C is a time average of the effective exhaust velocity, the relation between 
the fuel requirement of the total number of thrusts and the total velocity correc- 
tion is given by: 


L 

n= 1 


A v n 


= C Ln k, 


(14) 


where k denotes the mass ratio of the initial propellant mass to the final. The 
condition (14) holds for all mass ratios up to an exhaust velocity of one-tenth 
the velocity of light. 


ANALYTICAL TREATMENT 

The algebraic miss distance B at the target planet Jupiter is a function of 
time, position and velocity which remains constant between impulses. It is con- 
venient to suppose that the miss distance is positive if the spacecraft passes to 
the left of the target planet and negative if the spacecraft passes to the right of 
the target planet. 
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The control-effectiveness or the velocity-sensitivity vector BB/BV is to be 
evaluated along the nominal trajectory and its magnitude is generally a mono- 
tonic decreasing function of time. 

Before analyzing the sum » given by equation (13), let us make two as- 
sumptions about the main error distributions. 

The first assumption is that the distribution of thrust mechanization error 
is supposed to be normal with zero mean and standard deviation a. 

The second assumption is that the RMS velocity correction uncertainty in 
the direction of BB/ BV is negligible. 


Bearing in mind the above assumptions, the sum Sj^ can be written [3 ] as 
follows : 



B 0 = B(t 0 ) is the initial miss distance, which is a random variable, and sup- 
posed to be normally distributed with zero mean and standard deviation o- Bq . 


The first term of (15) corresponds to the case n = 1 and the second term 
is valid for n > 1. 

The miss distance I B | decreases after each correction and | b I - 0, on 
the average, as t - T , the arrival time. In other words, to an upper bound on 
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the expected values of the last | B n | corresponds an upper bound on the time to 
target at the last correction, so that 


> t* , 

where t x denotes the earliest time for the final correction. 


(17) 


Since |BB/BV | decreases to zero, the inequality (17) leads us to the following 
inequality 


B X _ TJ/fXX 

BV ‘ BV ’ B " B(t } 


(18) 


which gives an upper bound on BB,^ /BV. 
Minimization of S , or 

N 

2, . . .,N} occurs [3] when 


Minimization of S N , or the optimum spacing of correction times {t n ; n = 1, 


and that if 



BB X 


BV 

BV 


BB X 

bb 2 

^ b n-i 

BV 

BV 

BV 

bb 2 

” 3B 3 

3B » 

BV 

BV 

BV 


then t x =0, while if 




> 


= P 


BB ( 


a Bo >pa 


0 

BY 


< POr 


BV 


(19) 


( 20 ) 


( 21 ) 
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then t is obtained by 


3V pa 


( 22 ) 


Condition t = 0 states that if the launching error is sufficiently large in 
comparison with subsequent errors, the first correction to the spacecraft tra- 
jectory should be made in the vicinity of the earth, i.e., several days after launch 
so that one can get adequate early determination of miss distance. 

The criterion for choosing BB^/BV, which appears in (19), is [ 3 ]: 



Where R denotes the nearest approach to the Jupiter’s center, Figure 3, V R is 
velocity of approach relative to the Jupiter, and fx is Gauss’ constant for Jupiter. 
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Figure 3— Trajectory Near Jupiter (Hyperbola) 


10 


The total number of corrective thrusts N is given [ 3 ] by 


I 



(nearest integer), if a > pa 



if ct b, 


< pa 



The minimized F is 


(24) 



Actually the number of midcourse maneuvers in a time-open mission de- 
pends on many factors such as the characteristics of the nominal trajectory and 
objectives of mission, and so on. It can be seen from equation (24) that for a 
given nominal trajectory and an RMS launching error cr B ^ , the total number of 
midcourse maneuver N , though limited, increases with t N . 


NUMERICAL ANALYSIS AND RESULTS 

The nominal trajectory used for this galactic probe is adopted from [ 7 ] 
and has the following elements and characteristics: 

Heliocentric Eccentricity e = 0.8655 

Heliocentric Semi -major axis a = 1.09318 x 18 9 km 
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Heliocentric Semi-minor axis 
Heliocentric Semi-Latus rectum 
Launch date 
Injection time 
Geocentric injection speed 
Flight time to Jupiter 


b = 0.54758 x 10 9 km 
i = 0.27429 x 10« km 
December 30, 1969 
7 h 5 m 37 s GMT 
15.405 km/sec 
500 days 

M = 0.132715115 x 10 12 kmVsec 2 


The trajectory near the planet Jupiter (hyperbola) has the following 
characteristics: 


Eccentricity 

Semitransverse 

Radius of closest approach 
to Jupiter 

Velocity of approach rela- 
tive to the planet Jupiter 

^jupiter — ^* m jupiter 


= 2.479 

a t = -0.7077 x 10 6 km 

R = 10 6 km 

V R = 20.8 km/sec 

fi = 0.12671059 x 10 9 km 3 /sec 2 . 


In order to verify either of the inequalities (20), (21), and solve the equation 
(24), the numerical values of the standard deviations ct b , a, as well as the value 
of /3V, for a specified t , must be determined. 0 

Taking into account of the following numerical data [ 7 ] : 


cr ^ 2.4 x 10 6 km 
B o 

3B n 

~ 0.518 x 10 s km/m/sec 

BV - 

cr = 0 5x 10 -3 km/sec, 

one finds that the inequality (20) is satisfied, and therefore the first form for N 
in equation (24) should be adopted. 
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In order to get an estimate for the upper bound of t N one may refer to the 
guidance chart for Earth-Jupiter transfer [8] in Figure 4, which represents the 
variation of the control-effectiveness |BB /3V | along the above mentioned nomi- 
nal trajectory. 
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Figure 4-Guidance Chart for Earth-Jupiter Transfer 


From Figure 4 it can be seen that the approximate time for the upper bound 
of t N is 62 days from injection, and if we take for 

t N = 50 days 

after injection the optimum spacing ratio is 

P = 2.5 (26) 
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On the other hand 


= 5500 km/m/sec, 
BV 


therefore the optimum integer N , given by the first form of equation (24) is 

N = 3. (27) 

The sequence of times (t n > for which the condition (17) is satisfied will be 

tj = 10 days, t 2 = 35 days, t 3 = 50 days. (28) 

Figure 4 shows that | BB(t)/BV | is a monotonic decreasing function of time 
and the sequence 

BB. BB_ BB, 

= 35000 km/m/sec., = 14400 km/m/sec. , = 5500 km/m/sec. (29) 

BV BV BV 

is geometric with the approximate ratio 1/2.50. 

A glance at the literature shows that the optimum spacing ratio varies 
practically in the range 


2.1 < p < 3 (30) 

so that the value obtained for p is not substantially different from those obtained 
by others. 

The determination of fuel requirement for the total velocity correction as 
well as the magnitude of each velocity correction is of the utmost interest. In 
fact, the fuel expenditure is connected with the minimization of S N or F which 
is given by upper form of equation (25): 



where 

(1 + P 2 ) U2 / Ln p = 2.93857 for p = 2 . 50 . 
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As it can be seen from the above condition the first velocity correction AV t , 
depends on the first term in Eq. (31) , while the remaining corrections depend on 
logarithmic term. A detailed investigation of this problem should follow later. 


CONCLUSIONS 

The application of the spacing-ratio method to the case of a fly-by jovian 
mission leads us to the following conclusions: 

1. The optimum choice for the spacing of corrections is one which mini- 
mizes the total velocity correction or, equivalently, minimizes the 
total required fuel. 

2. A constant spacing-ratio of p * 2.5 for the velocity-sensitivity leads 
to a near optimum spacing of correction times. 

3. The total number of corrective thrusts obtained by this method is 
higher, and therefore it is advisable to use this method if higher accu- 
racies and more specific mission objectives, such as entry to the • 
atmosphere of a target planet, or soft landing, are required for the 
terminal trajectory. 
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